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Abstract 

The generalized Chaplygin gas (GCG) model explains the recent acceler- 
ated expansion of the Universe via an exotic background fluid whose equation 
of state is given by p = —A/p a , where A is a positive constant and < a < 1. 
The model is an interesting alternative to scenarios involving scalar field po- 
tentials, with the ensuing unnatural fine tuning conditions for the underlying 
particle physics theories. We derive constraints on the parameter space of 
the model from bounds on the location of the first few peaks and troughs of 
the the Cosmic Microwave Background Radiation (CMBR) power spectrum 
arising from recent WMAP and BOOMERanG data. 



PACS number(s): 98.80. Cq 

1 Introduction 

It has recently been proposed that the evidence for a dark energy component 
to the total energy density of the Universe at present might be explained 
by a change in the equation of state of the background fluid rather than 
by a cosmological constant or the dynamics of a scalar field rolling down a 
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potential PQ. This allows, at least in principle, to avoid well known fine- 
tuning problems associated with ACDM and quintessence models. Within 
the framework of Friedmann- Robert son- Walker cosmology, one considers an 
exotic background fluid, the GCG, which is described by the following equa- 
tion of state 



where a is a constant in the range < a < 1 (the Chaplygin gas corresponds 
to the case a = 1) and A a positive constant. Inserting this equation of state 
into the relativistic energy conservation equation, leads to a density evolving 
as [2] 



where a is the scale-factor of the Universe and B an integration constant. It 
is remarkable that the model interpolates between a universe dominated by 
dust and a De Sitter one via an intermediate phase which is a mixture of a 
cosmological constant and a perfect fluid with a "soft" matter equation of 
state, p = ap (a ^ 1) [7j. Notice that even though Eq. (JTJ) admits a wider 
range of positive a values, the chosen range ensures that the sound velocity 
(Cg = a A/ plh") does not exceed, in the "soft" equation of state phase, the 
velocity of light. Furthermore, as pointed out in Ref. J2|, it is only for 
< a < 1 that the analysis of the evolution of energy density fluctuations is 
meaningful. 

Furthermore, as discussed in Ref. j2], the model can be described by 
a complex scalar field whose action can be written as a generalized Born- 
Infeld action. Recently, it has been shown that a curvature self-interaction 
of the cosmic gas can mimic the GCG equation of state [3,. It is quite 
clear that the GCG is a candidate for explaining the observed accelerated 
expansion of the Universe [2j as it automatically leads to an asymptotic 
phase where the equation of state is dominated by a cosmological constant, 
ZttGA 1 / 1 ^. It has also been shown that the model admits, under conditions, 
an inhomogeneous generalization which can be regarded as a unification of 
dark matter and dark energy [21 E] without conflict with standard structure 
formation scenarios [21 El El- Hence, it is fair to conclude that the GCG model 




(2) 
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is an interesting alternative to models where the accelerated expansion of the 
Universe is explained via an uncancelled cosmological constant (see [7] and 
references therein) or a scalar field potential as in quintessence models with 
one |H] or two scalar fields [9j. Recently, some questions have been raised 
concerning the viability of the GCG model. For instance, in Ref. [TO], it 
is claimed that the model produces a matter power spectrum inconsistent 
with observation; however, the authors did not include the effect of baryons, 
which should play a crucial role and, in particular, would require a two-fluid 
analysis, as was done in Ref. [TT| . with the conclusion that the GCG can 
be quite different from the ACDM model and still reproduce 2dF large scale 
strucure data. On the other hand, in Ref. [TS], it is argued that the GCG 
model is indistinguishible from the ACDM model, which is not surprising 
as the authors did not consider the GCG as an entangled mixture of dark 
matter and dark energy as expected in a unification model. 

The possibility of describing dark energy via the GCG model has led to 
a wave of interest aiming to constrain the model using observational data, 
particularly those arising from SNe la O D3] and gravitational lensing 
statistics [To] . 

In this work, we extend the analysis carried out in Ref. (see also 
Ref. JH] for a study based on the CMBfast code) aiming to constrain the 
parameters of the GCG model from recent bounds on the positions of peaks 
and troughs of the CMBR power spectrum, employing basically the same 
methods that have been used to constrain quintessence models (see e.g. 
Refs. H2H2D1EUI22I)- Restricting the analysis of the CMBR power spectrum 
to the locations of peaks and troughs rather than considering the structure 
of the whole spectrum turns out to be a simple but very powerful tool in 
constraining the model parameters basically because of the precision with 
which these positions are now determined, especially following WMAP re- 
sults. We find that the model is compatible with WMAP bounds on the 
locations of the first two peaks and first trough, and BOOMERanG bounds 
on the location of the third peak provided a < 0.6, thus ruling out the 
Chaplygin gas model. The allowed range of model parameters depends, in 
particular, on h and n s ; for instance, for h = 0.71 and n s = 1, we obtain 
a < 0.4, 0.76 < A s < 0.88. These bounds become tighter for n s < 1 e.g. for 
n s = 0.93, we get a < 0.2, 0.79 < A s < 0.82. The allowed regions of model 
parameters become slightly larger for smaller values of h. 

Finally, we should like to mention that, in order to make the Chaplygin 
gas model consistent with the location of peaks and troughs in the CMBR 
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power spectrum as measured by WMAP, values of h smaller than the ones 
suggested by WMAP data are required, namely h < 0.65, together with the 
condition that n s is close to 1. 



2 CMBR constraints for the GCG model 

The CMBR peaks arise from acoustic oscillations of the primeval plasma just 
before the Universe becomes transparent. The angular momentum scale of 
the oscillations is set by the acoustic scale, Ia, which for a flat Universe is 
given by 

l A = 7T— , (3) 

c s r ls 

where r = J a _1 <it is the conformal time, r and T\ s being its value today 
and at last scattering respectively, while c s is the average sound speed before 
decoupling: 



Cs = r ls I c s dr , (4) 

where 



o 



with pb/p-y the ratio of baryon to photon energy density. 

In an idealised model of the primeval plasma, there is a simple relationship 
between the location of the m-th peak and the acoustic scale, namely l m ~ 
ml a- However, the peaks position is shifted by several effects which can be 
estimated by parametrising the location of the m-th peak, l m , as in [THJ 123 

V = Ia (m ~ Pm) = ^A(m -<p — 5ip m ) , (6) 

where (p = ipi is the overall peak shift and 5ip m = (p m — <p is the relative shift 
of the m-th peak relative to the first one. Eq. © can also be used for the 
position of troughs if one sets, m = 3/2 for the first trough and m = 5/2 
for the second trough. Even though analytical relationships between the 
cosmological parameters and the peak shifts are not available, one can use 
fitting formulae describing their dependence on these parameters. We use the 
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formulae given in Ref. j20] for the first three peaks and first trough, which 
we reproduce in the Appendix, for convenience. It is relevant pointing out 
that although these formulae were obtained for quintessence models with an 
exponential potential, they are expected to be fairly independent of the form 
of the potential and the nature of the late time acceleration mechanism as the 
shifts are practically independent of post recombination physics. We should 
stress that the analytic estimators we are using, determined by comparison 
with CMBfast for standard models, is less than one percent [20J. 

Following our dark matter-energy unification scenario, we rewrite the 
energy density, Eq. (J2J), as 



where A s = Aj p^wf anc ^ Pcho — {A + B) l l l+a . In terms of the new variables, 
Friedmann equation reads 



where we have included the contribution of radiation and baryons as these 
are not accounted for by the GCG equation of state. 

Several important features of Eq. (j7| are worth remarking. Firstly, we 
mention that A s must lie in the interval < A s < 1 as otherwise p c h would 
be undefined for some value of the scale-factor. Secondly, for A s = 0, the 
Chaplygin gas behaves as dust and, for A s — 1, it behaves like a cosmo logical 
constant. Notice that the Chaplygin gas corresponds to a ACDM model only 
for a = . Hence, for the chosen range of a, the GCG model is clearly 
different from the ACDM model. Another relevant issue is that the sound 
velocity of the fluid is given, at present, by aA s and thus aA s < 1. Using the 
fact that Pro/Pcho = £l r0 /(l - tt r0 - ^bo) and pw/Pcho = ^W(l _ ^ro ~ ^w), 
we obtain 







H 2 = n ch0 H*a~ 4 X 2 (a) , 



(9) 



with 



X(a) 




1 — fi r Q — fifeo 1 — ^rO — ^60 
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+ a 4 A s + 



a 3(l+a) 



1/1+as 



(10) 



Moreover, since H 2 = a 4 , we get 



dr 



da 



nHlH X(a) 



(11) 



so that 
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1 da I f a is da 



o X(a) \Jo X(a 



-i 



1 



'12) 



where ai s is the scale factor at last scattering, for which we use the fitting 
formula 



af, 1 -l = z ls = 1048 [1 + 0.00124V 738 ] [1 + <7i<] 



(13) 



where 



92 = 0.56[1 + 21.1^ 181 ]- 1 , (14) 

and UJb,m = Qb,mh 2 - 

Let us now turn to the discussion of the available CMBR data. The 
bounds on the locations of the first two peaks and the first trough, from 
WMAP measurements of the CMBR temperature angular power spectrum 
are 

£ P1 = 220.1 ±0.8, 
£ P2 = 546 ± 10 , 

£ dl = 411. 7 ±3.5; (15) 

where all uncertainties are la and include calibration and beam errors. The 
location of the third peak, from BOOMERanG measurements, is given by 

m 

i v , = 825+^0 . (16) 
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From the computation of the acoustic scale, Eq. (fT2j). the equation for the 
shift of the peaks, Eq. and the fitting formulae given in the Appendix, we 
look for the combination of GCG model parameters that is consistent with 
observational bounds. Our results are shown in Figures 1-4, where have 
assumed that uib = 0.0224 and used the fact that A s and Q m are related by 

1 — — f2 r 

which is obtained by noting that for a = the model is just the ACDM 
model; thus, one should identify the Chaplygin gas parameters with the 
usual density parameters when substituing a = in Eq. (8) (for the present, 
ao = 1), taking into account that Q m = Qb + ^cdm- 

In Fig. 1, we plot contours in the (h,Q m ) plane corresponding to the 
bounds on the first three peaks and first trough of the CMBR power spec- 
trum, Eqs. (fTljj) - (fTBj) . for n s = 0.97 and different values of a. The box 
on the a = (ACDM model) plot corresponds to the bounds on h and u m 
arising from the combination of WMAP data with other CMB experiments 
(ACBAR and CBI), 2dFGRS measurements and Lyman a forest data [23] : 

= 0.7118S , u m = 0.135j£°£ . (18) 

Notice that the above bound on h is slightly more restrictive than the bound 
obtained from WMAP data alone 



h = 0.72 ±0.05 . (19) 

Figures 2 and 3 show the same contours but for n s — 1 and 1.03, respec- 
tively. In Figs. 4 and 5, contours are shown in the (A s , a) plane for h = 0.71 
and h = 0.6. 



3 Discussion and Conclusions 

In this work we have shown that current bounds on the location of the first 
few peaks and troughs in the CMBR power spectrum, as determined from 
WMAP and BOOMERanG data, allow constraining a sizeable portion of 
the parameter space of the GCG model. Our results indicate that WMAP 
bounds imply that the Chaplygin gas model (a = 1 case) is ruled out and 
so are models with a > 0.6. For low values of n s , a = 0.6 is also ruled out. 
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However, for n s > 0.97, a = 0.6 becomes increasingly compatible with data. 
Hence, one can safely state that models with a < 0.2 are always consistent. 

Our analysis shows that results depend strongly on the Hubble parameter 
and since WMAP's bound on this quantity was obtained for ACDM models 
and, on the other hand, there are recent determinations of the Hubble con- 
stant, combining Sunayev-Zeldovich and X-ray flux measurements of clusters 
of galaxies, that give much lower values of H , namely |2Zj 

H = 60 ± 4+H km/s/Mpc , (20) 

and |2B| 

H = ffi±\* ± 15 km/s/Mpc , (21) 

it is relevant to examine the implications, in particular regarding the ex- 
clusion of the Chaplygin gas model, of relaxing the bound (|T9"j) and allow 
for lower values of h. Figures 4 and 5 show that, for h = 0.71 (the central 
value for WMAP's bound on h), a — 1 is not allowed for any combination of 
parameters; however, for h = 0.6 (slightly below WMAP's preferred range), 
a = 1 is allowed provided n s is around 1. In fact, a deeper analysis shows 
that, in order for the Chaplygin gas model to become consistent with peak 
and dip locations of the CMBR power spectrum, it is necessary that h < 0.65 
and n,wl, 

These results are compatible with the ones found in Ref. [T7] using bounds 
on the third peak from BOMERanG and the first peak from Archeops [29 
data as well as bounds from SNe la and distant quasar sources, namely 
0.2 < a < 0.6 and 0.81 < A s < 0.85. We find, in particular, that bounds 
from SNe la data, which suggest that 0.6 < A s < 0.85 [Hj, are consistent 
with our present results for n s = 1 and h = 0.71, namely 0.78 < A s < 0.87. 

Note added After we had completed this work, a related study has ap- 
peared jHOl which makes a likelihood analysis based on the full WMAP CMB 
data set using a modified CMBfast code, with results similar to ours. 
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Appendix 



We reproduce here the analytic approximations for the phase shifts found 
in Ref. [213 • The overall phase shift is given by 



where 



1.466 - OAQQn, 



a x r a * + 0.2910^ 



(22) 



ai = 0.286 + 0.626w 6 

a 2 = 0.1786 - 6.308^ + 174.9^ -1168^ (23) 
are fitting coefficients, 

% = < ^ 1 rn ch (r)dT, (24) 
J o 

and 

r*= Prad(zi s )/ 'pm{zis) (25) 

is the ratio of radiation to matter at decoupling and z\ s is given by Eqs. ()13j) 

-dH. 

There is no relative shift of the first acoustic peak, bxp\ = 0, and the 
relative shifts for the second and third peaks are given by 



Sif 2 = c - Cl r* - c 2 r; C3 + 0.05(n s - 1) , (26) 



where 



x 



r (1 = -0.1 + (0.213 - 0.123fif 

exp{- (52 -63.6fiJ?) uj b ) 



and 



ci = 0.015 + 0.063 exp (-3500^) , 
c 2 = 6 x Hr 6 + 0.137(c; b -0.07) 2 , 

c 3 = 0.8 + 2.30^ + (70 - 126^) u b , (27) 



5^3 = 10 - d x rf 2 + 0.08(n s - 1) , (28) 
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with 

d t = 9.97 + (3.3 - 3ft£) u b , 

d 2 = 0.0016 - 0.0067^ s h + (0.196 - 0.22^f) cu b 

+ (2.25 + 2.770$) x lO" 5 ^" 1 . (29) 

The relative shift of the first trough is given by 

5y?3/2 = b + hrl /3 exp(fc 2 r*) + 0.158(n s - 1) (30) 

with 

6 = -0.086 - 0.079^ - (2.22 - 18.10^) cu b 

- (140 + 403f&) ^ , 
= 0.39 - 0.98ft^ - (18.1 - 29.2ft£) u b 

+UO00I (31) 
b 2 = -0.57 - 3.8 exp(-2365c^ 2 ) , (32) 
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Figure 1: Contour plots of the first three Doppler peaks and first trough 
locations in the (fl m , h) plane for GCG model, with n s = 0.97, for different 
values of a. Full, dashed, dot-dashed and dotted contours correspond to 
observational bounds on, respectively, £ Pl , £ P2 , £ P3 and i^, see Eqs. (fTo^) and 
(I16|) . The box on the a = plot (corresponds to ACDM model) indicates the 
bounds on h and fl m h 2 from a combination of WMAP and other experiments, 
Eq. (JIEJ). 
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Figure 4: Contour plots of the first three Doppler peaks and first trough 
locations in the (A s , a) plane for GCG model, with h = 0.71, for different 
values of n s . Full, dashed, dot-dashed and dotted contours correspond to 
observational bounds on, respectively, £ P1 , £ P2 , £ P3 and i^, see Eqs. (fTo"]) and 

HI. 
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Figure 5: As for Fig. 4 but with h = 0.6. 
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